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2.  3 


Cylindrical  coordinates.  Fig.  1 

Longitudinal,  tangential  and  radial  displacements  of  the  cylindrical 
shell.  Fig.  1.  Note  that  a positive  displacement  w is  inward. 

Coefficients  defining  the  shape  of  modes  in  vacuo,  Eq.  (1) 

u,  v,  and  w components  of  displacement  pertaining  to  the  generalized 
coordinate  qk,  Eqs.  (2)  and  (3). 

Subscript 

Radius  of  cylinder 

Coefficient 

Velocity  of  sound  in  medium 

Element  of  surface  area  of  shell 

r -dependent  function  in  expression  for  potential 


yCI 

Length  of  longitudinal  half-wave  of  mode  of  vibration 
Generalized  mass  (coordinate  qk) 

Mass  of  shell  per  unit  area 
Virtual  mass  of  entrained  medium 
uw  » Nvw  Coefficients,  see  Eqs.  (40),  (46) 

Number  of  circumferential  waves  of  mode  of  vibration 
P(t,  9,  z)  - cos  n 9 sin  ^ e External  radial  force  (positive  if  outward) 

Pressure  in  medium 
Radial  pressure  of  medium  on  cylinder 
Static  pressure  in  medium 
Generalized  force  (coordinate  qk) 

Generalized  coordinate 
Large  value  of  the  radius  r 

Time  dependent  function  in  expression  for  potential 
Time 

w due  to  static  force 


Wk 


Uk2  + V*2  + WK2 


Coefficient 


- in- 


p,  (S  see  Eqs.  (19),  (37) 

A Coefficient,  see  Eqs.  (26),  (30),  (38)  and  Tables  1 and  2. 

p Mass  density  of  medium 

$=  S(t)F(r)  cos  n 0 sin  ^ Potential  function 

Q Frequency  of  free  or  forced  vibration  of  shell  surrounded  by  medium 

gjk  Frequency  of  cylinder  in  vacuo 
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Summary 


A method  is  presented  permitting  the  determination  of  the  frequencies  of  vibra- 
tions  of  infinitely  long  thin  cylindrical  shells  in  an  acoustic  medium.  (Sections  1 and  3) 

Expressions  are  also  obtained  for  the  displacements  of  the  shell  and  for  the  pres- 
sures p in  the  medium  in  the  case  of  forced  vibrations  due  to  sinusoidally  distributed 
radial  forces.  The  results  indicate  that  there  is  a low  frequency  range  where  no  radi- 
ation takes  place,  and  a high  frequency  range  where  the  external  force  provides  energy 
which  is  radiated.  Resonance  occurs  in  the  low  frequencies  range  only,  in  the  high 
range  it  is  prevented  by  the  damping  due  to  radiation. 

Free  and  forced  vibrations  of  steel  shells  submerged  in  water  are  discussed  in 
detail  (Sections  3 and  4).  With  limitations,  the  theory  may  be  applied  approximately  to 
stiffened  shells. 

The  method  requires  only  a minor  modification  to  account  for  the  effect  of  static 
pressure  in  the  surrounding  medium  (Section  5). 

The  results  are  also  used  to  draw  a general  conclusion  concerning  the  treatment 
of  transient  problems,  (Section  6),  particularly  of  the  problem  of  shock  loads- 
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i.  Free  Vibrations 


The  free  and  forced  vibrations  of  a submerged  infinitely  long,  thin  cylindrical 
shell  are  studied  in  this  paper  by  considering  the  shell  without  the  fluid  as  a separate 
structure  responding  to  the  applied  forces  combined  with  the  dynamic  forces  exerted  by 
the  surrounding  infinite  acoustic  medium.  Using  the  modes  of  free  vibration  of  the  shell 
in  vacuo  as  generalized  coordinates,  its  response  can  be  expressed  in  terms  of  the  in- 
finite number  of  these  modes. 

The  frequencies  and  modes  of  vibrations  of  infinitely  long  thin  cylindrical  shells 
required  for  the  purpose  were  determined  in  ref.  1.  The  individual  modes  can  be  clas- 
sified by  the  length  L of  the  longitudinal  half-wave  of  the  displacements  and  by  the  inte- 
gral number  n of  circumferential  waves.  For  each  length  L and  number  n there  exist 
three  frequencies  o)k(k  = 1,  2,  3).  The  displacements  of  the  modes  corresponding  to  all 
three  frequencies  are  sinusoidal,  but  the  ratios  of  longitudinal,  circumferential  and 
radial  displacements,  u,  v and  w differ.  Excluding  the  case  n = 0,  the  displacements 
corresponding  to  the  frequency  wk  may  be  written: 

u = Uk  cos  n 9 cos  ^5. 

v = Vk  sin  n 9 sin  ^ (1) 


w = Wk  cos  n 9 sin  -j- 


The  modes  defined  in  Eqs.  (1)  contain  an  arbitrary  factor,  only  the  ratios  Vk  /Uk  and 
Wk/Uk  can  be  determined  and  are  contained  in  ref.  1.  It  might  also  be  mentioned  that 
Eqs.  (1)  can  be  generalized  by  adding  phase  angles,  which  is  equivalent  to  changing  the 
origin  of  the  coordinates  9 and  z. 

The  most  general  response  of  a cylindrical  shell  will  be  a combination  of  the 


modes  for  all  values  of  L and  n.  However,  in  order  to  find  for  a submerged  shell  the 
frequencies  and  modes  of  free  vibration  of  selected  values  L and  n,  one  need  consider 
only  the  three  modes  of  shape  (1)  having  the  same  parameters  L and  n.  This  simplifi- 
cation is  due  to  the  fact  that  the  wave  equation  (8)  for  the  surrounding  medium  has  a 

7TZ 

solution  containing  the  factor  cos  n 9 sin  -j—  , and  the  pressure  belonging  to  such  a so- 
lution docs  not  excite  modes  of  different  values  L or  n. 

Selecting  the  coefficients  Wk  in  Eqs.  (1)  as  generalized  coordinates  Wk  = qk,  and 
using  the  abbreviations 


Uk 

= w; 

Vk 

Wi; 


cos  n 9 cos  -y- 

. n . TTZ 

sin  n 0 sin  -j- 


(2) 


cos  n 0 sin  -j- 


the  displacements  u,  v and  w can  be  expressed  as  functions  of  the  three  generalized  co- 
ordinates q| , q2,  q3 


u = qi  4>u,i 

+ qa4>u,2 

+ q3<t>ut3 

(3a) 

V = qi  4>v,i 

+ q24>v,2 

+ qs^v.s 

(3b) 

w = (q,  + q2  + q3H>* 

(3c) 

The  symbol  4>w  in  Eqs.  (2)  and  (3)  does  not  have  a second  subscript,  as  there  is  no  dif- 
ference in  the  w - displacements  of  the  three  coordinates  qK. 

Due  to  the  inherent  orthogonality  between  the  displacements  represented  by  the 
generalized  coordinates,  the  equations  of  motion  contain  each  only  one  of  the  coordinates 


qk. 


M„qk  + Mkiokqk  = Qk 


(k  = 1,  2,  3) 


(4) 


where  the  generalized  mass  Mk  is 


Mk 


+ 4>J' 


9w  )dA 


(5) 


In  the  case  of  free  vibrations  the  generalized  force  Qk  is  due  solely  to  the  radial  pres- 
sure pa  of  the  medium  on  the  cylinder  at  r = a: 


Hence,  Eq.  (4)  becomes 


Qk  = 


Po.<K»dA 


(°) 


qK 


2 

wk  3k 


11 


Pa<MA 


m (4>u,k  + 4>v,k  +4>«)dA 


(7) 


The  integrals  in  Eq.  (5)  to  (7)  are  to  be  taken  over  the  surface  of  a section  of  the  cylin- 
drical shell  of  length  L. 

The  shell  is  assumed  to  be  submerged  *n  -.a  acoustic  medium,  the  velocity  poten- 
tial $ for  which  is  governed  by  the  equation 


V3*  .-L  «± 

c2  at2 


(8) 


A solution  of  this  equation  suitable  for  the  present  problem  is 


$ = S(t)F(r)  cos  n € sin  ^ = S(t)F(r)<|>w 


(9) 


where  3(t)  and  F(r)  are  functions  at  our  disposal.  At  the  surface  of  the  shell  (r  = a),  the 
radial  velocity  of  the  shell  must  equal  that  oi  the  medium: 


-2- 


w = 


8$ 

dr  J r„a 


Substituting  Eqs , (3c)  and  (9)  one  obtains 

(q  i + q2  + q3)<t>w  = - S(t)F'  ( a)<t>w 


(10) 


or 

- (q,  + q2  + q3) 

S<t)  prjjj (11) 

(note:  F' (a)  - SfiEl]  _ J. 

The  radial  pressure  exerted  on  the  shell  by  the  surrounding  medium  is  the  pres- 
sure at  r = a: 


Pa  * " P §Tjr.a  = " P^C^)^r(a) ‘t’w 


(12) 


or  substituting  Eq.  (11), 


Eqs.  (6)  and  (7)  become 


Pa  = (qi  + P2  + Qs) 


(13) 


Qk  = p + ^2 + qs)  Jj^wdA 


(14) 


and 


where 


2 Pak  Ffal 

qk  + “kqk  = -jjp  pf^a)  (qi  + q2  + <b) 


(15) 


11 


4>wZdA 


WiT 


o-k 


||  (4>u,lf  + <t>v,kZ+  <t>w)dA  Uk  + Vr  T 'Vk 


(16) 


The  last  equation  is  obtained  by  substituting  Eq.  (2).  Numerical  values  of  a.k  for  various 
modes  are  given  in  the  tables  of  ref.  1. 

The  function  F(r)  is  obtained  from  Eq.  (8),  which  reads  in  full 


a2$  j_  a$  j_  §j|_  8^$_  _ j_  e^_ 

~~T  Z 2 + . 2 2 2 

9r  r ar  r 30  9z  c at 

Assuming  periodic  solutions  for  the  generalized  coordinates  qk 

„ int 
qk  = Cke 


(8a) 


(17) 


Eqs,  (9)  and  11)  give 


Introducing  the  symbol 


* = ' inFi75vt*»»  (Ci  + c?.  + ^s)e 


L c 

Eqs.  (8a)  and  (18)  lead  to  the  differential  equation 


F"(r)  + p F'(r)  -(p2  + I^F(r)  = 0 


As  boundary  conditions  we  require  that  the  velocity  of  the  medium  at  r = » vanishes  in 
such  a manner  that  the  kinetic  energy  in  the  medium  remains  finite. 


A solution  of  Eq.  (20)  is 


r"  F(r)  L.oe  " 0 


F(r)  = H‘”  (ipr) 


which  satisfies  Eq.  (21)  provided  the  real  part  of  p is  positive. 
Computing 

F'(r)  = ip[ip-HBC°  (ipr)  + HB.,l°(ipr)] 


Eq.  (15)  becomes: 


where 


F(a)  . 

Tfa-T" 


n - ipa 


Hf-V’  (ipa) 
Hn”  (ipa) 


3k  + “k^k  _ “ ak  A(3|  + ^2  + <b) 


* = *<Pa>  = X F® 


H n-i  (ipa) 

n - ipa  _Tn 

H*  (ipa) 


Substituting  Eqs.  (17)  into  (25)  three  homogeneous  equations  in  the  coefficients  C k are 
obtained 

(«k2- n2)Ck  - (C,  +C2  + C3)  = 0 (k  = 1,  2,  3)  (2' 


-4- 


f 


Non-vanishing  solutions  of  these  equations  exist  only  if  their  determinant  vanishes, 
which  after  evaluation  leads  to  the  frequency  equation 


o>i2  - n2 


a2 

cj22-  QZ 


“3  _ 2m 

2 2 2 
wj  - ^ paQ  A 


(28) 


valid  for  n f 0 . 

For  the  previously  excluded  case,  n = 0,  the  cylindrical  shell  in  vacuo  has  again 
three  modes  for  each  value  L and  n.  One  of  these  modes  is,  however,  a purely  tor- 
sional motion  with  displacement  components  u = w = 0.  This  mode  is  not  excited  by 
the  radial  forces  from  the  surrounding  medium,  and  only  the  two  non -torsional  modes 
need  be  considered.  Instead  of  Eq.  (25)  one  obtains  for  n = 0; 


Qk  + Ulfok  = "■^^■nkMqi  + ^2) 


(29) 


where 


iH9tl!  <i(3a) 

A = - 2 

pa  H|  (ipa) 

The  frequency  equation  for  n - 0 is; 


(30) 


az 


- + 


2m 


2 _ 2 2 _2  — .2 

c*)  1 - S2  - Q pa  Q A 


(31) 


The  solution  of  the  frequency  equations  (28)  and  (31)  is  further  discussed  in  Section 


3. 


2,  Forced  Vibrations  due  to  Radial  Forces 

Any  arbitrarily  distributed  radial  force  P(f,  0,  z)  can  be  expanded  in  a Fourier 
integral  (or  series)  of  terms 

P(t,  0,  z)  = ein*  cos  n 0 sin  - eln,<j>w  (32) 

The  following  treatment  will  be  restricted  to  the  consideration  of  such  individual  terms, 
i.e.  to  force  components  P(t,  0,  z)  defined  by  Eq.  (32).  These  forces  P are  counted 
positive  if  acting  outwards.  Non -radial  forces  can  be  treated  in  a similar  manner  but 
will  not  be  considered. 

Eqs.  (4)  and  (5)  of  Section  1 apply  again,  but  Eq.  (6)  for  the  generalized  force  Qk 
must  be  replaced  by 


: «k 


ini  2 
e 4>w  dA 


(33) 
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Proceeding  essentially  in  the  same  manner  as  in  Section  using  Eqs.  (8)  to  (13) 
incl„  and  Eq,  (16)  to  (24)  inel.,  the  equations  of  motion  for  n f 0 become 


qh  +a,l!2qh  = -|^aKA(q,  +q2  + q3)  -^ein<  (k  = 1.  Z,  3)  (34) 

There  is  just  one  detail  in  the  derivation  of  this  equation  which  requires  elucida- 
tion. The  solution  (22) 

F(r)  = HB‘"(ipr) 

TTC 

does  not  satisfy  the  boundary  condition  (21)  for  frequencies  n > -jj-  because  the  argument 

ap'a^  '71 

becomes  purely  imaginary.  The  boundary  condition  (21)  is  in  this  case  replaced  by  the 
requirement  that  the  potential  function  Eq.  (18) 

* = 'FTSy(Ci  +C*  +c3)eint  FW-  (35) 

represents  an  outgoing  wave  for  large  values  of  r.  This  condition  is  satisfied  for 

F(r)  • h[2)  (pr)  (36) 

where 


and  the  positive  sign  of  the  square  root  applies.  This  leads  to  the  value 

2 


A = A(?a)  = * §T^)  = 


n - f)a  • 


Hn.|2'  (Pa) 
Hn‘2’  (fa) 


(38) 


valid  for  while  for  Q4  Eq,  (26)  remains  in  force. 

Substitution  of  Eq.  (17) 


q*  = Ck  e 


int 


in  the  equations  of  motion  (34),  leads  to  three  non-homogeneous  linear  equations  for 
the  coefficient  Ck.  Their  solution  is 


Ct, 


a« 

2 o2 

wk  - n 


(39) 


where 


-6- 


(40) 


Ql  $7 


Qo  fl 


N = — 


a3f2‘ 


2 _2  2 „2  2 

ui  - n t»2  - n cjo g - q 


Substituting  the  values  Ck  from  Eq.  (39)  into  the  expression  (18)  for  the  potential 
function  $ and  noting  Eqs.  (26),  (32)  and  (3B)  leads  to 


i flaA  » 
2m 


L2m 


o2  a n 

mn  A ' ir_ 


e.  «) 


(41) 


where  the  function  F(r)  is  given  by  Eqs.  (22)  and  (36),  respectively, 


c H„(l)  (ipr) 


if  n < 


TTC 


F(r)  *< 


if  ^ and  n f 0 


^ H„!2)  (fTr) 


if  n > 


TTC 


(42a) 
(42b) 
(47  cl 


A,  0 and  5 are  defined  by  Eqs.  (19),  (26),  (37)  and  (38), 
The  pressure  in  the  medium,  p - - p * becomes 

1 F(r) 


p = 1 __  2rn  F(a) 


P(t,  0,  z) 


(43) 


The  pressure  pa  on  the  surface  of  the  cylinder  is 

Ptt  = 


i - 


dm 

paTSTA 


(43a) 


The  knowledge  of  the  coefficients  Ck  permits  also  the  determination  of  the  dis- 
placements of  the  cylinder.  Of  particular  interest  are  the  radial  displacements  w 
obtained  from  Eq,  (3c) 


...  . 1 P(t,  9,  z) 


m n 


Pa  * . 

ZSA  TT 


(44) 


If  one  goes  to  the  limit  n-*0  this  equation  gives  an  expression  for  the  deflection 
w*t  under  static  loading  P(9,  z)  = 4>w 


wit  = lim  w = - 
Q — o 


U|  Qg  Q3" 

+ + 

2 2 2 
■ CO  j CO  2 CO3  - 


m 


(44a) 


while  a similar  process  for  the  pressure  leads  to  the  obvious  result  lim  p = 0, 

O-o 

The  expressions  for  the  longitudinal  and  tangential  displacements  u and  v due  to 
the  radial  force  P(t,  0,  z)  are  somewhat  more  complicated 
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Ti 


Nuw 

u - 

cos  n 0 cos  ^e10* 

2 

mN  Q 

pa  A . 1 
1ST*  TT 

N 

nvw 

V = 

sin  n 6 sin  lnt 

mN  ft2 

Pa  A - 1 

2 N 

where 


SL-J  2 T 

k = 1 Wk  - » 


^ 


1^,  Qk  wkn 


/ J 2 2 

ifTi  u*  - ft 

The  above  Eqs.  (40)  to  (46)  were  derived  for  the  case  n f 0,  but  with  exception  of 
Eq  (42b),  they  remain  valid  for  n = 0 if  the  terms  referring  to  the  non -existing  frequency 
u3  are  omitted  when  computing  N,  Nuw,  and  N»„  from  Eqs.  (40)  to  (46). 

The  degenerate  solutions  (42b)  for  n = 0 and  for  n - 1 do  not  satisfy  Eq.  (21)  and 
the  total  kinetic  energy  in  the  medium  is  in  these  cases  not  finite.  This  indicates  that 
steady  state  vibrations  for  n = 0 or  1 in  an  actually  infinite  medium  are  not  possible  if 
ft  = . If  one  considers  the  infinite  boundary  of  the  medium  as  an  approximation  of  a 

far  boundary,  r = R,  the  solution  (42b)  for  n = 1 may  be  considered  an  approximation  of 
physical  meaning  because  the  energy  is  then  finite  and  the  velocities  at  the  boundary 
are  small,  as  1/R.  A similar  reasoning  does  not  apply  for  n - 0,  because  the  longitudinal 
velocity  is  not  small  for  large  values  of  r;  the  physical  impossibility  of  this  solu- 
tion expresses  itself  in  the  fact  that  A - °°.  To  obtain  a physical  meaningful  solution  for 
n = 0 it  is  necessary  to  consider  the  actual  conditions  at  the  distant  boundary. 

To  facilitate  numerical  computations,  Tables  1 and  2 list  the  values  of  A occuring 
in  Eqs.  (41)  to  (46)  for  n = 0 to  4 incl.* 


3,  Discussion  of  the  Frequency  Equation 
The  discussion  can  be  restricted  to  the  frequency  equation 


I'Values  of  A for  n ^ 20  can  be  obtained  from  ref.  4,  but  only  if  ft  > • The  acoustic- 

reactance  and  resistance  ratios_Xn  and  0„  shown  in  Fig.  1 and  2 of  ref.  4 can  be  used 
to  computem^A^)=  Xn(ap)  - i9„(ap). 
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and  examination  sbj  „ ■ , , 

„c  ;>ws  that  the  -values  ot  t\  are  real,  it  aj3  is  real,  but  that  a becomes 
complex,  if  $"2  > -» — , 

*-  so  lh&;  ».p  iu  imaginary:.  This  le ;;dp  to  the  question  how  many  real 
hivnish,  Che  coefficients  ths  frequencies  «k  on  the  left  hand 

2tl! 

pa 

’ solution  i Vi  whur,  ap  ,1s  imaginary,  that  is  any  real  root 
Q must  be  below  the  lin’d  ••  ‘ 


roots  ft  Eq.  (20)  wil 
side  of  Eq.  (23)  being  , 
real  Eq,  (28)  can  have  * 


e o \ ,‘he  <»•■•  tire  left  hand  side  will  always  be  real.  As  is  also 


Further  information 


Eq.  (28)  graphically.  For  n , 


n the  number  of  real  roots  o.n  be  obtained  by  considering 


of  the  left  hand  side  oi  c,q.  (2i 


< 0,  Fig.  2 shows  in  solid  lines  a typical  plot  of  the  values 


' as  function  of  ft.  Starring  with  a positive  value  for 


< 


iymptot.es  for  each  of  the  vaiues  ii  * wi ,vz  and  wj  . The 
►,  ted  ja  dashed  lines,  ti  e intersections  indicating  roots 


TTC 


1 right  hand  side  has  real  valuer,  only  for  0 4 ft  .<  -j~ 
lot  indicates  that  thery  are  at  least  as  many  real  roots 
the  system  in  vacuo  which  lie  below  -■■■■ , and  at  *bf 
> can  occur  if  sheer  ia  a further  frequency  w*  j»'jt  tar*«j 


ft  = 0,  the  plot  shows  vertical  a 
right  hand  side  of  Eq.  (28)  is  ' 
of*he  frequency  equation.  As 
the 'curve  ends  at  j oii»r  E.  The  ? 
ft  as  there  are  frequencies  w*  of 
most  one  more.  T:  is  latter  cas, 

than  so  that  ths  endpoint  frt’  . ... 

E Is  in  r..  position  lido  E1  '.h  Fig,  t. . u in  parcvinUr  id, 

frequencies  wi,  are  '-a  rtf  .r  than  ■ .c  ....  , 

• there  sny  be  no  lea-.  ireq^ency  at  ail.  ir-c  case  oi 
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For  n = 0 somewhat  sin/ 

/.iar  conditions  prevail,  Fip.  >.  *n  this  case  the  *.-.••«  /wu 
Eii.es  however  on  thw  ‘.7  - ax?  .....  r - , . ,, 

,,.  I which  results  in  the  occurence  of  ».ne  real  r.;vt  jn  >U 

P VP'fl  if  Lj . Is--  ^ - . A 27  / 

* ’ •■*.«  ;xa.mp\e  of  such  a case  is  a oteeJ  shell  in  water  tnrraliiiff 

with  a half  wave  ien.qtn  I,  > ij 


It  is  also  a*>p  .treat  fi  •' 


tion  mu£f  be  small  h-  Uian  j 


im  Figs  2 and  3 that  the  lowed  root  r>(  of  the  frequency  eq* 


between  u,  ond  ui  and  (? 
above  the  highest  freque'i* 
Attentio  n is  drawn  . 

sncU  in  vacuo 


.he  frequency  w,  , while  the  next  root  fij,  , if  any,  mu  at  lie 
ie  highest  root  must  be  between  and  «3  . No  roof  can  lie 
j u>j  in  vscuo, 

■)  the  fact  that  due  Ic  the  orthogonality  of  the  modes  q,  of  the 


U( 


Q2  + (13  = 1 


(4ft) 
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As  a rule  one,  and  for  many  configurations  two  values  a*  are  quite  small,  indicating 
that  the  radial  displacements  of  such  modes  are  small.  In  such  oases  the  coupling  of 

;7  I f’  i; 

the  motion  of  the  shell  and  the  surrounding  medium  is  slight,  and  the  frequencies  £1  of 
such  modes  will  be  very  close  to  the  frequencies  u in  vacuo. 

To  simplify  the  numerical  determination  of  the  roots  of  the  frequency  equation 
Table  1,  giving  values  of  a for  real  arguments  ap,  has  been  provided.  After  a frequency 
rt  has  been  obtained  the  shape  of  the  mode  can  be  determined  by  computing  the  coefficient 
C i,  from  Eqs,  (27).  ’ 

Effect  of  Compressibility 

It  is  of  nuereut  to  consider  the  effect  of  the  compressibility  inforder 

to  determine  when  it  might  be  neglected.  The  compressibility  appears  in  the  derivation 
of  Eq.  (8)  where  the  velocity  of  Refund  c occurs,  and  the  solution  for  an  incompressible 
fluid  may  be  obtained  by  using  c . This  does  not  change  the  frequency  equatioq  (28), 
but  the  argument  ap  used  in  deter  (tuning  the  term  A in  this  equation  is  in  this  casje 


an  - av 

aP  T 


(49) 


instead  of  the  value  given  in  Eq.  (47)  which  depended  on  the  frequency  fi  . Comparing 
Eqs.  (47)  and  (49),  it  is  apparent  tfhat  the  compressibility  may  be  neglected  only  if 


Q-«i 

c U 


If  the  compressibility  is  neglepted  the  number  of  real  roots  Q is  always  equal  to 


the  number  of  frequencies  <ak . it  U‘<  obvious  that  roots  higher  than  ^ found  for  th  i in 

( * .01  ° . ij 

compressible  case  are  unreal.  Mi 


Complex  Roots  of  Frequency  Equation  1 

The  fact  that  the  number  of  rear,  roots  n of  the  frequency  equation  may  be  /Tier 
than  the  number  of  frequencies  w*  in  vacuo  raises  the  question  of  the  existence  oi  ( qm- 
plex  roots  in  lieu  of  the  missing  real  foots,  and  of  the  physical  meaning  of  comafet 
roots  if  any.  In  order  to  have  physic  ,1  meaning  the  solution  corresponding  to  a ?o*|i  '.  V3X 
root  n must  represent  an  outgoing  wa-  e,  and  mist  s’ so  decay  with  time. 

: c.:. 

In  order  to  recognize  the  wave  character  at  the  solution  it  is  convenient  to 
potential  function  in  the  form  (36)  employed  in  flection  2 for  Q > , 

Ce;mHna)  (j5r)*w  (50) 

To  avoid  incoming  waves  the  real  pa*-;  • >f  must  not  be  negative.  Based  on  this  expres- 
sion one  obtains  a "wave  frequency  equation"  . which  is  identical  with  Eqs.  (?Ji)  or  (31), 
except  that  the  value  A is  a function  of  ap  anc5  must  be  computed  from  Eq.  (38)  instead 
of  Eq.  (26).  It  is  to  be  noted  that  the  «rgumsn‘:  aP  in  Eq.  (38)  will  now  be  complex.  A 
study  of  the  potential  function  (50)  indicates  that  for  complex  roots  Q for  which  the  time 
dependent  term  e mt  is  damped,  the  to.  ;.  H„ (pr)  becomes  oo  at  r -»  °o  . Such  roots  and 
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the  corresponding  solutions  appear  at  first  to  have  no  physical  significance.  It  should 
be  pointed  out  however,  that  these  complex  roots  and  the  corresponding  solutions  occur 
and  have  significance  as  asymptotic  solutions  in  transient  problems.  In  such  cases  the 
validity  of  the  asymptotic  solution  is  limited  to  finite  values  of  r and  the  fact  that 
H„  ((?r)  does  become  infinite  at  r -»oo  is  immaterial.  For  example,  a heavy  shell  in  a 
tight  nr|ed  um,  like  air,  is  capable  of  slightly  damped  vibrations.  The  frequency  equa- 
tion in  such  a case  has  complex  roots  f2  the  real  parts  of  which  are  close  to  the  frequen- 
cies ii  vacuo,  while  the  imaginary  parts  are  quite  small.  Expression  (50)  is  then 
the  approximation  of  a possible  motion,  valid  in  the  range  t > T and  r < R,  where  T 
must  be  selected  so  large  that  the  transient  effects  of  the  initiation  of  the  motion  have 
passed  ou.side  the  cylinder  having  the  arbitrarily  selected  radius  R. 

The  numor*  aal  values  vf  complex  frequencies  n can  easily  be  determined  approxi- 
mately either'  the  imaginary  part  of  n is  small,  or  if  101  is  large.  If  the  imaginary 
part  of  fl  s sma”  the  value  of  A(a{J)  computed  for  the  adjoining  real  part  of  Q,  found  in 
Table  2.  may  be.used.  This  procedure  would  be  applicable  for  the  above  mentioned  case 
of  a heavy  shell  In  a light  medium,  if  the  absolute  value  of  the  frequency  |fll  is  large 
such  that  the  ar,  ament  lajJI  » 1,  the  term  A(ap)  may  be  determined  by  using  asymptotic 
expressic  is  for  the  Hankel  Functions, 


A(af)  = j + 2.ap 


(51) 


Frequently  it  u;  permissible  to  neglect  even  the  real  part  of  the  denominator  in  this 
equation,  .jiving. 

2i 


= - 


a? 


(52) 


These  apy.  -oxirv 

of  n.  n ix.  •*«*»'•*  c'  i>Y. 
suits  of  f v Si/  w 

Virtual 


>?.,U  values  of  n,  but  give  better  results  for  small  values 
1 ecJ.  of  a”'  'oximation  can  be  obtained  by  comparing  the  re- 
••':h  vah  ci  for  real  arguments  af£  in  Table  2. 


Fur  j ■.  :u  ..e  ooss.h.Hs  to  simplify  the  frequency  equation  by  neglecting  the  two 
higher  thus  obtuw/g  an  approximate  equation  for  the  lowest  frequency.  This 

equation  can  luievyre.ed  in  r :<i?.u'ple  manner  using  the  concept  of  the  " virtual 

mass"  of  hit*  «•»’  t valued  medium. 

values  of  an  and  uk  occuring  in  Eq.  (6b;  are  giv«p  in  ref.  1.  These 
,Ai..  0 : ‘"'•yd thin  the  range  1^  10  the  first  term  of  Eq.  (28),  a(/(W|  - O*  ),  As 


'•'tjther  two  if  n < wi  j for  — values  near  1 this  is  so  because 
1 “g  2 

?nd  a3  a-..:  ^ina’!  fur  i3v>>  *lues  because  «2  and  ws  are  very  much  larger  than 


very  much  ?r./ver  that 

id  ay 

Negle'-iivf,  the  iwns  w.Uh  £ “*! 


d <*3  leads  after  rearrangement  to 


a,  pa 


i:  (1  + 15T  A)  = 


(53) 
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This  indicates  that  the  lowest  frequency  n of  the  submerged  shell  is  equal  to  the  fre- 
quency of  a shell  in  vacuo  whose  mass  has  been  increased  by  the  " virtual  mass"  m»  of 
a portion  of  the  medium  vibrating  with  it 

m,  = a,  ^ A (54) 

A is  a function  of  the  frequency  given  by  Eqs.  (19)  and  (26),  but  if  wi  « the  effect  of 
the  frequency  is  negligible  and  Eq.  (49)  for  the  incompressible  case  may  be  used  instead 
of  Eq.  (19).  mv  may  be  computed  quickly  by  the  use  of  Table  1. 

The  same  approach  can  be  used  also  for  n = 0,  but  gives  good  results  only  if  — < 2. 
For  longer  shells  the  second  term  in  Eq.  (51)  remains  large  and  can  not  be  neglected, 
thus  making  the  simple  virtual  mass  concept  inapplicable. 

Steel  Shells  Submerged  in  Water 


When  applying  the  frequency  equations  to  thin  steel  shells  (say  r > 30)  submerged 

L 1 

in  an  infinite  body  of  water  it  is  found  that  for  ratios  1 < — 4 10  and  n > 1 one  and  only 

ire  a 

one  frequency  in  vacuo  is  lower  than  -j-  , while  for  n = 0 all  frequencies  wk  are  above 
this  limit.*  It  follows  that  for  n = 0 one  and  only  one  real  root  £2  exists;  for  n > 0 there 
is  at  least  one  root  and  according  to  the  previous  discussion  the  possibility  of  a second 
one.  A second  root  n can  only  occur,  see  E'  in  Fig.  3,  if  the  second  frequency  «2  is 
close  to  ^ . Ref.  1 shows  that  in  the  range  1 < ^ ^10  the  second  frequency  «2  is  much 
larger  than  the  limit  ^ , and  there  will  therefore  be  just  one  real  frequency  for  each 


larger  than  the  limit  , and  there  will  therefore  be  just  one  real  frequency  for  each 
value  of  n.**  This  frequency  can  be  determined  from  the  simplified  Eq.  (53),  except  if 
n = 0 and  — > 2. 

cl 

Similar  conditions  prevail  for  ratios  L/a  > 10,  except  for  n » 2,  where  for 
extremely  large  ratios,  due  to  the  effect  of  the  bending  stiffness  of  the  shell  all  the 
frequencies  in  vacuo  may  lie  above  ^ , resulting  in  no  real  roots.  In  the  limiting 

case  L = <»  only  complex  roots  occur  for  any  value  of  n.  This  case  has  been  discussed 

•! 

independently  in  ref.  2. 

If  — < 1 the  above  result  is  expected  to  remain  valid,  unless  the  length  of  the  half 

:<  & 

wave  is  so  short,  comparable  to  the  thickness  t,  that  all  frequencies  in  vacuo  lie  above 

TT  • 

Complex  roots  of  the  frequency  equation  representing  waves  exist  and  can  be 
determined  approximately  by  using  the  asymptotic  values  for  A(aj5),  Eq.  (52).  However 
the  solutions  corresponding  to  the  complex  roots  lack  physical  meaning,  except  as 
asymptotic  solutions  of  transient  problems  discussed  in  the  previous  section.  In  the 


♦ The  frequencies  in  vacuo  for  1 ^ < 10  can  be  obtained  from  ref.  1. 

♦♦This  count  excludes  the  frequency  of  the  purely  torsional  mode  of  the  shell  for  n = 0. 
This  mode  was  excluded  when  deriving  Eq.  (31)  as  it  is  not  coupled  with  the  surround- 
ing medium. 
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case  of  a steel  shell  in  water  where  1 < ^ 4 10,  the  asymptotic  solution  is  a combination 
of  the  free  vibration,  corresponding  to  the  real  frequency  , and  of  damped  wave  mo- 
tions corresponding  to  the  complex  roots  Q . If  the  motion  is  such  that  the  undamped 
motion  is  at  all  excited,  the  asymptotic  solution  converges  rapidly  toward  the  undamped 
vibration  of  real  frequency  fl  and  the  damped  motions  due  to  complex  frequencies  will 
give  only  unimportant  contributions, 

4.  Discussion  of  Results  for  Forced  Vibrations 
Response  of  Shell  and  Medium 

The  most  significant  feature  of  the  results  of  Section  2,  is  the  fact  that  the  re- 
sponse is  entirely  different  for  frequencies  fl  of  the  applied  force  above  and  below  the 
limit  fl  = . Different  expressions,  see  Eqs.  (42  a - c),  apply  in  the  two  ranges,  re- 

sulting in  different  physical  behavior. 

If  f2<  ^ the  ratios  Hr,111  (i(3r)/Hnli'  (ipa),  and  as  a consequence  the  values  of  the 
coefficient  A(a|3),  are  real.  Eqs.  (43)  to  (45)  indicate  that  in  this  case  the  pressure  p 
and  the  displacements  are  in  phase  with  the  applied  force  P,  There  is  no  input  of 
energy  by  the  applied  force  and  no  radiation.  It  is  of  interest  that  the  pressure  p de- 
creases very  rapidly  with  increasing  radius  because  HI!')  (ipr)  decreases  like 
e~pr  jyr.  As  expected  the  expressions  for  the  pressure  and  displacements  contain  a 
denominator  which  vanishes  if  becomes  equal  to  one  of  the  roots  of  the  frequency 
equation  for  free  vibrations,  indicating  resonance 

If  the  frequency  Q is  sufficiently  smaller  than  -j-  the  effect  of  the  compressibility 
becomes  unimportant  and  the  approximate  Eq.  (49)  for  ap  may  be  used.  For  low  fre- 
quencies, to,  one  can  neglect  the  modes  qz  and  qj  . The  resulting  equation  can  be 
interpreted  by  means  of  the  virtual  mass  defined  by  Eq.  (54),  indicating  that  the  response 
is  identical  with  that  of  a shell  of  increased  mass  (m  + m, ) in  vacuo. 

If  Q > conditions  are  quite  different.  The  ratio  HnZ)  (pr)/felnl2)  (pa)  is  complex 
and  the  coefficient  A(ap)  also  is  complex.  As  a result,  the  pressure  p and  the  displace- 
ments are  not  in  phase  with  the  applied  force,  the  phase  angle  being  sucn  that  the  ap- 
plied force  supplies  energy  which  is  radiated.  There  are  no  real  roots  of  the  frequency 
equation  above  ^ , and  resonance  does  not  occur  in  this  range. 

The  range  Q > can  be  further  subdivided  in  two  distinct  ranges  if  n / 0.  Table 
2 indicates  that  for  small  values  of  ap  the  imaginary  part  of  A is  very  much  smaller  than 
the  real  one,  while  for  large  values  of  aP  the  opposite  is  true.  The  dividing  line  is  quite 
sharp,  but  depends  on  the  value  of  n.  For  n 2.  2 it  lies  above  aP  = n/2.  This  indicates 
that,  particularly  if  n >2,  there  is  an  extensive  range  of  frequencies  ft  > ^ where  ra- 
diation and  damping  are  quite  small.  The  decay  of  the  pressure  p as  function  of  r is 
also  different  in  the  two  ranges.  If  the  radiation  is  large  the  pressure  decays  very 
slowly,  approximately  as  ]/  yr.  If  the  radiation  is  small,  the  pressure  decays  much 
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faster,  nearly  as  fast  as  in  the  non-radiating  case  (although  for  very  large  values  of  r 
the  decay  becomes  again  1/  yT). 

Considering  the  response  as  function  of  the  forcing  frequency  n.  a number  of 
significant  frequencies,  in  addition  to  frequencies  of  free  vibrations,  can  be  found  by 
considering  the  denominators  of  Eqs.  (43)  to  (45).  If  the  frequency  O is  equal  to  one  of 
the  frequencies  uk  in  vacuo,  N = » , and  Eq.  (44)  gives  pa  = P(t,  9,  z);  the  pressure  in 
the  medium  is  in  this  case  equal  to  the  applied  pressure,  as  if  the  shell  were  non- 
existent. This  is  physically  understandable,  because  at  such  a frequency  the  shell 
vibrates  without  requiring  any  outside  force,  and  the  eritire  applied  force  is  transmitted 
to  the  medium. 

Other  significant  values  of  the  frequency  n are  those  for  which  the  term  N vanishes. 
Consideration  of  Eqs.  (40)  and  (46)  shows  that  this  occurs,  in  addition  to  the  trivial 
value  0 = 0,  once  between  successive  values  of  the  frequencies  wk  in  vacuo.  If  N = 0, 
both  p and  w vanish  everywhere,  indicating  that  there  is  no  response  at  all  of  the  fluid 
at  such  frequencies.  This  rather  startling  result  is  due  to  the  fact  that  for  these  fre- 
quencies the  three  modes  (or  two  modes  if  n = 0)  interact  in  the  manner  of  a vibration 
damper.  If  the  longitudinal  and  radial  displacements  u and  v are  determined  it  is  seen 
that  they  do  not  vanish,  only  w and  p are  zero.  The  frequencies  at  which  this  phenome- 
non occurs  are  necessarily  the  same  as  for  the  shell  in  vacuo. 

At  the  end  of  Section  2 it  has  been  pointed  out  that,  and  for  what  reason,  the  solu- 
tion for  n = 0 loses  its  meaning  for  n = ^ . The  formulas  for  the  response  for  frequen- 
cies very  close  to  n = ^ must  also  be  used  with  caution,  because  they  might  be  affected 
by  the  conditions  at  the  distant  boundaries  of  the  medium. 

Steel  Shells  in  Water 

Considering  cases  in  the  range  1 4 ^ 10,  there  is  for  n 1 one  frequency  in 

vacuo,  w,  , below  ^ , while  the  other  two,  w2,  w3.  are  much  higher.  The  frequency 
equation  of  the  submerged  shell  has  one  real  root,  Oi  < and  two  complex  ones  near 
and  <i)3  . 

Figs.  4 and  5 show  the  amplitudes  of  the  radial  response  w and  the  maximum 

pressure  pa  at  the  surface  of  the  shell  in  non-dimensional  manner  for  a typical  case, 

n = 2,  t = 4,  JiiL  s 4,  (The  ratio  ^ = 4 indicates  that  the  mass  of  the  water  displaced 
am'  m 

by  the  shell  is  2.0  times  the  mass  of  the  shell)  Fig.  4 shows  in  solid  lines  the  ratio* 

j w/w stat f c j of  the  amplitude  of  the  displacement  w to  the  static  deflection  of  the  shell 

as  function  of  the  ratio-^1  . Attention  is  drawn  to  the  fact  that  the  right  half  of  the 

ire 

diagram  is  enlarged  and  shows  100|w/to,to,|C|  . The  same  figure  shows  for  comparison 
in  dashed  lines  the  response  of  the  shell  in  vacuo.  Fig.  5 shows  the  ratio  of  the  peak 
values  of  the  surface  pressure  pa  to  the  peak  value  P of  the  applied  pressure  P(t,  0,  s). 


♦Often  called  amplification. 
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The  diagrams  show  clearly  the  previously  mentioned  fact  that  w and  p become 

zero  at  two  points,  where  N = 0,  which  is  due  to  the  vibration  damper  effect.  The  dif- 

or 

ference  between  the  response  in  vacuo  and  in  water  for  ^ 1 as  shown  in  Fig.  4 is 
essentially  that  between  an  undamped  and  damped  system,  the  response  curve  showing 
peaks  near  the  undamped  natural  frequencies  wi,  . 


The  principal  point  of  interest  in  the  response  curve  for  the  pressure. 


Pa 

T" 


Fig. 


5,  is  its  sawtooth  pattern.  There  is  the  expected  resonance  at  the  one  real  root 


followed  by  a minor  minimum' near  = 1;  then,  near  <*>2  and  W3,  there  are  zero  values 
followed  closely  by  peaks  of  unit  value.  It  is  also  of  interest  that  the  pressure  decreases 
with  increasing  frequencies  only  slowly,  much  slower  than  the  response  w, 

Figs.  6 and  7 show  the  radial  response  w and  the  pressure  p4  for  the  case  n = 0. 

The  difference  between  the  responses  w in  vacuo  and  in  water  is  much  more  pronounced 
as  in  the  case  n = 2.  This  is  principally  due  to  the  fact  that  the  only  frequency  V.  i of  the 
submerged  shell  is  very  much  lower  than  the  fundamental  frequency  u|  of  the  shell  in 
vacuo.  Attention  is  drawn  to  the  extreme  narrowness  of  the  peaks  of  the  curves  in  Figs. 
6 and  7 at  the  resonant  frequency  fi ! . The  response  curve  w near  = i has  been 
omitted,  because  at  this  point  the  solution  for  an  infinite  medium  becomes  meaningless 
(See  the  last  paragraph  but  one  of  Section  2). 


Stiffened  Shells 

It  is  worth  noting  that  the  theory  developed  is  applicable,  approximately,  to  ring 
stiffened  shells  provided  the  length  L of  the  half  wave  is  several  times  the  stiffener 
spacing.  In  such  a case  modes  of  vibrations  of  the  stiffened  shell  exist  for  which  the 
displacements  are  approximately  sinusoidal,  so  that  the  modes  may  be  expressed  again 
in  the  form  of  Eq.  (1).  If  the  frequencies  and  the  shapes  of  the  modes  in  vacuo  are 
known*  the  previously  derived  formulae  can  be  applied  and  furnish  approximate  results 
for  stiffened  shells.  Such  a treatment  should  give  good  approximations  at  low  frequen- 
cies, while  at  high  ones  detail  will  be  lost,  because  the  vibrations  of  the  shell  between 
stiffeners  have  been  neglected. 


5.  Effect  of  Static  Pressure 

The  theory  presented  in  Sections  1 and  2 neglected  the  effect  of  a static  pressure 
in  the  surrounding  medium  on  the  shell. 

The  effect  of  the  static  pressure  on  the  medium  is  indirectly  allowed  for  by  the 
use  of  an  appropriate  value  c for  the  velocity  of  sound,  but  the  "buckling  effect"  of  the 
static  pressure  on  the  shell  does  not  appear  in  the  previous  derivations. 

Consider  the  shell  in  vacuo  on  which  an  external  pressure  p0  (smaller  than  the 


*It  is  intended  to  treat  engineering  methods  for  the  determination  of  these  modes  in  a 
future  report. 
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buckling  one)  is  applied  by  an  imaginary  inertialess  device.  The  modes  of  free  vibra- 
tions and  frequencies  can  be  determined,  furnishing  a set  of  new  modes  having  in 
general  lower  frequencies  than  in  the  case  p0  = 0.  To  account  for  the  effect  of  the 
static  pressure  it  is  necessary  to  use  these  new  modes  and  their  frequencies  in  the 
previous  analysis. 

Provided  p0  is  smaller  (say  at  least  10%)  than  buckling  pressure  of  the  shell  in 
the  mode  having  the  same  half  wave  length  L and  circumferential  wave  number  n,  the 
shape  of  the  modes  in  vacuo  for  zeito  pressure  and  for  a pressure  p0  will  not  differ 
substantially.  One  can  conclude  that  this  is  the,  case,  by  applying  the  reasoning  used  ih 
ref.  1 (pp.  28,  29)  to  show  that  the  shape  of  the  modes  in  a thin  shell  will  not  be  appre- 
ciably affected  by  its  bending  stiffness.  It  is  therefore  permissible  to  use  the  shapes  of 
the  modes  without  static  pressure  (given  in  ref.  1)  as  approximation  of  the  new  modes, 
but  new,  corrected  values  iop  for  the  frequencies  must  be  used.  The  corrected  frequen- 
cies cop  can  be  found  from  Raleigh's  principle.  If  this  is  done  it  is  found  that  even  the 
frequencies  change  only  very  slightly,  except  the  fundamental  ones  for  n 1.  It  is 
therefore  sufficient  to  correct  the  fundamental  frequencies  for  n ^ 1,  and  use  the  origi- 
nal ones  for  n = 0,  and  for  the  higher  frequencies  if  n > 1. 

The  additional  term  in  the  expression  for  the  potential  energy  due  to  a radial 
pressure  p0  is  known  from  the  theory  of  buckling  (ref.  3,  p.  130) 


w 


2 2 t 

a - u « - 2awu,  Jdfidz 


(55) 


where  the  symbols  are  those  of  ref.  1.  To  be  able  to  apply  the  theory  for  the  infinitely 
long  shell  to  a shell  of  great,  but  finite  length  it  is  of  interest  to  be  able  to  include  also 
the  effect  of  a longitudinal  compression  due  to  the  pressure  p0  on  the  end  surfaces.  The 
additional  contribution  to  the  potential  energy  is 


Po 

T 


+ v. 


]d9dz 


(55a) 


Applying  Raleigh’s  principle  the  corrected  fundamental  frequency  is 


\f~2  Po 
“P=  ~amB 


where,  in  the  case  of  radial  pressure  only 


B = 


(n2  - 1)  W,2  - 2 W,  U | + n U,Z 
- 5 o 2 

U,  + V,  + w, 


(56) 


(57) 


while  for  combined  radial  and  longitudinal  pressures 


-16- 


, 1 \ W ^ 11/  TT  tt  2 1 tt  ^ , 1 3L  7T  2 

(n  - 1)  Wi  -2  -y—  W|  U|  + n U|  + •*■  V|  + T — 5— W| 

B = h r L_g_.  (57a) 

U|2+  V|2  + w(2 

The  values  of  B have  been  determined  for  n = 1 to  4 incl.  and  are  listed  in  Table  3 for 
1 < ^ 10.  It  will  be  noted  that  some  of  the  values  B for  n = 1 are  negative,  indicating 

an  increase  of  the  frequency  due  to  radial  pressure  only. 

6.  A Conclusion  Concerning  Transient  Problems 

The  results  of  the  analysis  of  . forced  vibrations  may  be  utilized  for  the  solution 
of  transient  problems.  In  such  cases  the  force  can  be  expanded  in  multiple  Fourier 
series  (or  integrals),  the  individual  terms  of  which  will  have  the  response  found  in  Sec- 
tion 2. 

The  fact  that  there  is  a fundamental  difference  in  the  response  depending  on 
whether  f2  is  larger  or  smaller  than  , may  be  used  to  predict  the  character  of  the 
response  of  transient  problems,  and  especially  the  type  of  approximation  permissible 
for  a particular  problem. 

In  shock  and  impact  problems  an  important  part  of  the  force  will  be  in  the  range 
of  high  frequencies.  The  response  during  and  shortly  after  the  application  of  the  force 
will  be  governed  by  the  compressibility  of  the  medium  and  radiation  will  be  important. 

If  on  the  other  hand  th~  response  is  to  be  studied  after  a certain  time  has  elapsed  the 
high  frequency  effects  will  have  vanished  and  the  elastic  and  inertial  effects  will  be  para- 
mount. 

Consider,  for  example,  the  case  of  the  transverse  bending  stresses  due  to  shock 
or  impact  in  a submerged  shell  or  other  structure.*  In  order  for  these  stresses  to 
reach  their  maximum,  a time  comparable  to  the  fundamental  period  of  vibration  must 
have  elapsed.  In  this  time  the  compressible  effects  had  time  to  decay,  and  the  elastic 
and  intertial  ones  control.  It  is  therefore  possible  to  study  the  response  with  good  ap- 
proximation by  an  expansion  into  the  modes  of  free  vibration  of  the  submerged  structure. 
It  should,  however,  be  remembered  that  the  results  obtained  in  such  a manner  are  not 
valid  for  the  pressure  in  the  medium  or  other  effects  shortly  after  the  impact. 

A case  where  the  compressible  effects  are  paramount  is  the  response  of  a cylin- 
drical shell  to  a transverse  step  shock  wave  treated  in  ref.  2.  The  major  effect  occurs 
in  the  two  lowest  modes  within  a time  much  too  short  for  the  compressible  effects  to 
decay.  In  this  case  the  response  of  the  shell  was  expressed  in  terms  of  the  modes  of 
the  structure  in  vacuo,  treating  all  the  fluid  pressures  as  external  forces  acting  on  the 
shell.  The  use  of  the  modes  in  vacuo,  and  not  of  the  modes  of  the  submerged  shell,  is 


*The  paper  considered  cylindrical  shells,  but  it  is  evident  that  the  general  behavior  of 
all  types  of  structures  must  be  similar. 


required  in  cases  involving  infinite  media  where  the  compressibility  remains  important. 
This  is  due  to  the  fact  that  for  a shell  in  an  infinite  medium  for  each  value  of  n only  a 
finite  number  of  frequencies  and  modes  exist,  even  if  complex  frequencies  are  included . 
(In  the  above  mentioned  example  there  is  for  n = 0 just  one  frequency  which  is  complex.) 
This  finite  number  of  modes  is  not  sufficient  to  express  the  infinite  number  of  states 
(for  each  value  of  n)  of  which  the  system  consisting  of  shell  and  medium  is  capable. 

The  number  of  modes  available  is  sufficient  to  express  any  state  of  displacement  and 
velocity  of  the  shell  at  a given  moment,  but  additional  terms  representing  waves  in  the 
medium  of  appropriate  pressure  and  velocity  must  be  superimposed  to  express  all  pos- 
sible simultaneous  states  of  the  medium.  The  use  of  the  modes  of  the  submerged  struc- 
ture is  in  such  a case  of  no  advantage  as  the  additional  terms  mentioned  remain  of 
paramount  importance.  If  in  a shock  problem  sufficient  time  has  elapsed  for  the  effects 
represented  by  these  terms  to  decay  the  solution  converges  asymptotically  towards  the 
modes.  Because  of  the  decay  of  the  modes  belonging  to  complex  frequencies,  only  the 
free  vibrations  of  real  frequency  will  remain.  This  is  of  course  nothing  else  but  the 
already  discussed  case  where  compressibility  is  unimportant. 

It  can  therefore  be  stated  that,  if  high  frequency  terms  occur  in  the  force,  or 
shock  effects  are  wanted  within  a short  time  after  the  application  of  the  force,  a treat- 
ment using  solely  modes  of  vibration  of  the  submerged  structure  would  be  incomplete, 
as  additional  terms  occur  in  the  solution.*  As  an  alternative  approach  the  modes  of 
free  vibration  of  the  structure  may  be  used  as  generalized  coordinates  describing  the 
response  of  the  structure  fully,  but  leaving  the  medium  to  be  treated  by  means  of  the 
differential  equations  for  the  potential,  or  in  any  other  way  desired. 
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Table  1 


Values  of  A if  Q < 


- % 
L c 

n = 0 

n = 1 

n = 2 

n = 3 

n = 4 

0 

00 

2.  000 

1. 000 

0.  667 

0.  500 

0.  1 

4.  928 

1. 952 

0.  997 

0.  666 

0.  4998 

o.  2 

3.  670 

1. 863 

0.  990 

0.  664 

0.  499 

0.  4 

2.  551 

1.  661 

0.  965 

0.  658 

0.  497 

0.  6 

1. 989 

1.473 

0 029 

0.  648 

0.  493 

0.  8 

1.  640 

1.  312 

0.  888 

0.  635 

0.  487 

1.0 

1. 399 

1.177 

0.  844 

0,  619 

0.481 

1.  5 

1.  028 

0.  928 

0.  737 

0.  575 

0.  460 

2.  0 

0.  814 

0.  761 

0.  645 

0.  528 

0.  436 

{ 3.0 

0.  577 

0.  556 

0.  505 

0.  443 

0.  385 

4.  0 

0.  447 

0.  437 

0.  411 

0.  375 

0.  338 

5.  0 

0.  365 

0.  360 

0.  344 

0.  322 

0.  298  : 

6.  0 

0.  309 

0.  305 

0.  296 

0.  281 

0.  264 

7.  0 

0.  267 

0.  265 

0.  259 

0.  249 

0.  237 

| 8.0 

0.  236 

0.  234 

0.  230 

0.  223 

0.  214 

9.  0 

0.  211 

0.  210 

0.206 

0.  201 

0.  195 

10.  0 

0.191 

0. 190 

0.187 

0.184 

0.  179 

Assymptotic  solution  for  all  values  of  n 


A 


4 

2 pa  + 1 


Table  2 
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Table  3 


Values  of  B in  Eq.  (56) 


Radial  Pressure  Only 


Combined  Radial  and  Longitudinal  Pressure 


1. 

00 

1. 

25 

1. 

50 

1. 

75 

2. 

00 

2. 

25 

2. 

50 

2. 

75 

3. 

00 

3. 

50 

4. 

00 

5. 

00 

6. 

00 

7. 

00 

8. 

00 

n 

= 1 

0. 

282 

0. 

111 

-0. 

0427 

-0. 

0457 

-0. 

212 

-0. 

240 

-0. 

248 

-0. 

245 

-0. 

236 

-0. 

211 

-0. 

185 

-0. 

140 

-0. 

107 

-0. 

0842 

-0. 

0672 

-0. 

0547 

2.  604 
2.362 
2.  222 
2. 152 
2.  126 

2.  123 
2.  132 
2.  147 
2.  165 
2.200 

2.  231 
2.  278 
2.  310 
2.  330 
2.  345 

2.  356 
2.  364 


7. 

108 

6. 

956 

6. 

906 

6. 

903 

6. 

922 

6. 

947 

6. 

973 

6. 

997 

7. 

019 

7. 

056 

7. 

083 

7. 

120 

7. 

143 

7. 

157 

7. 

166 

11. 718 
9.869 
8.923 
8.391 
8.071 

7.865 
7.  726 
7.  629 
7.559 
7.466 

7.410 
7.  348 
7.  316 
7.298 
7.286 
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